Hysteretic Depinning and Dynamical Melting for Magnetically Interacting Vortices in 

Disordered Layered Superconductors 
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We examine the depinning transitions and the temperature versus driving force phase diagram for 
magnetically interacting pancake vortices in layered superconductors. For strong disorder the initial 
depinning is plastic followed by a sharp hysteretic transition to a 3D ordered state for increasing 
driving force. Our results are in good agreement with theoretical predictions for driven anisotropic 
charge density wave systems. We also show that a temperature induced peak effect in the critical 
current occurs due to the onset of plasticity between the layers. 



A vortex lattice subject to quenched disorder offers an 
ideal system for studying glassy dynamics governed by 
the competition between the ordering forces (the vortex- 
vortex repulsion) and disorder jjj . One of the key issues 
of glassy dynamics is the nature of the depinning tran- 
sition separating a pinned state of the lattice at small 
external drives from the sliding regime that occurs when 
the external force exceeds a threshold depinning value. 
The early numerical simulations revealed that in the sys- 
tem with strong disorder plastic depinning occurs with 
the vortices tearing past one another and flowing in com- 
plex patterns . The onset of plasticity at depinning was 
later confirmed in experiment || as well as subsequent 
numerical studies and was demonstrated to play a 
crucial role in the pre-melting peak effect where the crit- 
ical current exhibits a striking increase as a function of 
temperature or field |l]|j| prior to vortex lattice melting 
and related phenomena. As the applied driving force is 
further increased, the effect of the quenched disorder is 
reduced and the vortex lattice can show a remarkable 
reordering transition at a driving force F r [Q. At this 
transition a considerable portion of the vortex lattice or- 
der is restored and the motion is highly coherent. The 
reordering force F r is expected to diverge as the temper- 
ature T approaches the vortex lattice melting temper- 
ature T m Q. Transport measurements provide strong 
evidence for this scenario ||]. Additional studies on the 
depinning and reordering transition and the nature of 
the coherently moving phases have been conducted the- 
oretically H , numerically Q] , and experimentally [p|^0| , 
showing both non-hysteretic continuous depinning transi- 
tions and abrupt hysteretic V(I) curves. The appearance 
of the hysteretic portion of the V(I) curves depends on 
how the vortex lattice was prepared: the initial current 
ramp shows a large hysteretic critical current but subse- 
quent ramps produce only continuous non-hysteretic V(I) 
characteristics Q. 

The underlying mechanisms of depinning and its inter- 
play with the reordering transition remain an open ques- 



tion. The issues to be clarified include a determination 
of the conditions under which the depinning transition 
shows hysteretic or switching responses, as well as an 
understanding of the nature of these phenomena. The 
mechanism by which a three-dimensional lattice depins 
or reorders is also an open question, as well as the issue 
of whether the onset of plasticity coincides with a sudden 
jump in the critical current (the peak effect). 

Some of these issues were addressed in pdf| in a context 
of a driven three-dimensional anisotropic charge density 
wave system (CDW) with disorder. It was shown that if 
disorder is sufficiently strong jl^] , the dependence of the 
lattice velocity on the applied drive becomes iwo-valued 
or bistable: at the same applied force (close to the de- 
pinning force of the topologically disordered phase) the 
lattice can either remain at rest or slide depending on 
the history. This means that the depinning transition 
shows switching behavior and hysteresis. The details of 
depinning are determined by particulars of the system 
involved. For a layered anisotropic CDW considered in 
|]llf the depinning occurs in two stages. The initial plastic 
depinning, where the decoupled 2D CDWs slide indepen- 
dently in each layer, is followed by a sharp transition or 
crossover at higher drives to a faster coherently moving 
3D solid CDW phase. Reversing the process reveals hys- 
teresis: pinning (immobilizing) of a coherent 3D CDW 
occurs at smaller drives. It was also shown that for weak 
disorder there is only a single non-hysteretic depinning 
transition directly to the 3D coherently moving state. 

The results of flUl are very general, and the particu- 
lar model used there can be straightforwardly general- 
ized to other periodic media in quenched disorder. A 
detailed study of the character of the many- valued v — F 
dependence near depinning was recently carried out in 
Jl3| within the framework of a mean-field model for a 
visco-elastically driven vortex configuration, where a dy- 
namic bistability and, accordingly, a hysteretic plastic 
depinning was found to occur for sufficiently strong dis- 
order. Note that the layered CDW model of [O can 
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be conveniently applied directly to highly anisotropic su- 
perconductors where the magnetic interactions between 
pancake vortices dominate (m|. In this system both a 
2D regime, in which vortices in each layer move inde- 
pendently from the other layers, and a 3D regime, where 
vortices in different layers form coherently moving lines, 
should exist 

In this work we investigate the dynamics of magneti- 
cally interacting vortices in a three-dimensional layered 
superconductor to numerically test the predictions of Rcf. 
jnj . Our model is relevant to anisotropic layered su- 
perconductors such as BSCCO as well as artificially lay- 
ered superconductors in which the magnetic interactions 
dominate. Recent numerical simulations with this model 
have found that for increasing applied driving force, a 
transition from a decoupled plastically flowing state to a 
coupled 3D coherently moving lattice occurs Jll|[l6|]. A 
sharp static 2D to 3D disorder induced transition also 
appears as a function of field and disorder [^5|JT^]. In 
this model, the dynamic phase diagram as a function of 
driving force and temperature, as well as hysteretic and 
switching behavior in the current- voltage characteristics, 
have not been studied until now. We find that for strong 
disorder, the initial depinning is 2D with the vortices 
flowing plastically and independently in any one layer. 
For increased drives we see a sharp transition to a co- 
herently moving 3D vortex lattice which manifests itself 
in an abrupt increase in the vortex velocities. We ob- 
serve strong hysteresis when the driving current is cycled. 
We have also investigated the driving versus temperature 
phase diagram. 

We consider a three-dimensional layered superconduc- 
tor in which the vortices are modeled as repulsive point 
particles confined to a layer with an equal number of vor- 
tices per layer. The overdamped equation of motion for 
vortex i is given by 

U = -£ ViU{pij, zv) + C + U + f T - VVi , 

where N v is the number of vortices, pij and z^j are the 
distance between pancakes i and j in cylindrical coordi- 
nates, and Tj = 1. The system has periodic boundaries 
in-plane and open boundaries in the z direction Q . The 
magnetic energy between pancakes is |H| , [l9f 

U( Pi ,,0) = 2 rf eo (V^)-f + ^) 



U(pij,z) = -s m — — f exp(-z/A) In — - E 2 j 

where eo = &q/ (4ir\) 2 , A is the London pen- 
etration depth, R = 22. 6A, the maximum ra- 
dial distance, E\ = dp' exp(— p' / A)/ p' , E2 — 



J p °° dp' exp(- y/z 2 + p' 2 /A) /p', d = 0.005A is the inter- 
layer spacing, and we set s m = 2.0. The uncorrelated 
pins are modeled by parabolic traps that are randomly 
distributed in each layer. The vortex-pin interaction is 

given by ff = -E£f(/p/&)fa - ^Mtiv - |r, - 
rj^D/A), where the pin radius £ p = 0.2 A, the pinning 
force f p = 0.02/q, and /q = eo/A. We have simulated a 
16A x 16A system with a vortex density of n v — 0.35/A 2 
and a pin density of n p = 1.0/ A 2 in each of L = 8 layers. 
This corresponds to N v = 80 vortices and N p = 256 pins 
per layer, with a total of 640 pancake vortices. For our 
finite temperature simulations the temperature is mod- 
eled as a Langevan kick f T where < f(i) T >= and 
<f(t)ft(t) T >=2r,k B T8 ij . 

In Fig. 1(a) we show the average vortex velocities in 
the direction of drive, V x —< v x >, as the applied driving 
force fd is increased and then decreased, and in Fig. 1(b) 
we show the corresponding c-axis correlation function, 
C z = 1 - ((|(r i)£ - r i)i+ i)|/(ao/2))e(ao/2 - \{r iiL - 
r i.L+i)|))j where eto is the vortex lattice constant. At 
fd = the vortices are pinned and uncorrelated in the 
z-direction as is indicated by the fact that C z has a low 
value. The initial depinning, at fd P = 0.01/q, is 2D with 
the vortices remaining uncorrelated in the z-direction. 
Near fd — 0.013 the vortex velocities sharply increase, 
coinciding with a sharp transition to a coupled moving 
3D state as seen in the jump up in C z p|. At higher- 
drives the vortex velocity linearly increases with increas- 
ing fd- Upon reducing fd from this linear regime the sys- 
tem remains in the 3D state for driving forces well below 
fd = 0.013. At fd = 0.007 there is a sudden transition to 
a 2D pinned state as seen in the drop in V x and C z . If the 
applied driving force is again increased the velocities and 
C z will follow the same hysteresis loop. In Fig. l(c,d) 
we show V x and C z for the same system as in Fig. l(a,b) 
but with a lower disorder strength of f p = 0.01/g. Here 
the depinning transition is non-hysteretic with a single 
depinning threshold to the 3D coherently moving state 
as seen by the fact that C z = 1.0 at all times. For a 
given coupling strength value s m , there is a critical pin- 
ning force above which there is a static 3D-2D transition 
Jl5t - Hysteretic V(I) curves can only be observed above, 
but close to, this transition. The width of the hysteresis 
in the V(I) curves is largest for disorder strengths just 
above the static 2D-3D transition and gradually narrows 
for increasing disorder strength. 

The behavior we observe is in excellent agreement with 
the theoretical predictions of We note that unlike 
the anisotropic charge-density wave system, where only 
inter-plane plasticity or slipping can occur, in the vortex 
system in-plane plasticity is possible within each layer. 
The 2D-3D transition in the case of the vortices occurs 
by simultaneous recoupling of the vortices between layers 
and the reordering of the vortices in plane. To illustrate 
this, in Fig. 2(a) we show a top view of the vortices in 
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FIG. 1. (a) Vortex velocities V x versus fd for a system with 
f p — 0.02/q at T = as the driving force is first increased 
from fd = 0.0 to fd = 0.02 and then decreased back to fd = 0. 
(b) C z versus fd for the system in (a) showing the sharp, hys- 
teretic recoupling and decoupling transitions, (c) V x versus 
fd for a system with f p = 0.01/o. (d) C z versus fd for the sys- 
tem in (c) showing that there is a non-hysteretic 3D pinned 
to 3D moving transition. 

the moving 2D regime, indicating that the vortices are 
uncorrelated from layer to layer and are disordered in 
plane. In contrast, Fig. 2(b) shows that in the moving 
3D regime, the vortices are aligned between layers and 
possess a triangular ordering. 

In Fig. 3 we show the driving force versus tempera- 
ture phase diagram for a system with f p = 0.01/q. The 
temperature is plotted in units of T m — 0.00045, the tem- 
perature at which the clean system undergoes melting in 
the form of a single sharp transition from a 3D lattice to a 
2D pancake gas. As shown in Fig. 3, at T/T m < 0.005 fol- 
low drives there is a 3D pinned phase where the vortices 
remain aligned. At these temperatures, for fd/ fo — 0.001 
there is a non-hysteretic depinning transition directly to 
the coherently moving 3D state. 

We find a static transition to the 2D state at T/T m = 
0.005, which we label Td c - For T > Td c the depinning 
transition is 2D, and as the driving force is further in- 
creased a dynamically induced reordering of the vortices 
occurs. As T is further increased above T^ c the size of 
the 2D pinned phase is reduced, while the drive at which 
the 2D plastic flow to 3D coherent flow transition oc- 
curs diverges as T approaches T mi in agreement with the 
dynamical freezing model of Koshelev and Vinokur Q . 

The thermally induced decoupling seen in Fig. 3 also 
coincides with a sharp increase in the critical current or 
a "peak effect." There has been considerable work on 
dimensional transitions in the context of the peak effect 
in layered superconductors. The static 3D-2D transition 
and critical current increase in layered superconductors 
have been previously studied as a function of interlayer 
coupling and magnetic field |2fJ,El| . Koshelev and Kes 
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FIG. 2. Pancake vortex positions seen from the top of the 
sample, for the system with f p = 0.02/o shown in Fig. l(a,b). 
Pancakes on a given layer are represented by circles with a 
fixed radius; the radii increase from the top layer to the bot- 
tom. Vortices that have aligned into vortex lines appear as 
thick circles, (a) At fd = O.Oll/o, in the 2D plastic flow 
regime, (b) At fd = 0.015/q in the coherently moving 3D 
phase. 

have proposed theoretically that magnetically coupled 
vortices can show a first order transition from a 3D to 
a 2D system under certain conditions p2[ |. The sharp- 
ness in the 3D-2D transition we observe in our simulation 
suggests that it is first order in nature p^ , p"7| |. 

The 3D-2D transition can be seen as a temperature 
induced peak effect in which the combination of temper- 
ature and quenched disorder bring about a decoupling 
transition where the individual pancakes can adjust to 
the optimal pinning configurations. The phase diagram 
in Fig. 3 is very similar to the one proposed by Bhat- 
tacharya and Higgins || in which there is an elastic de- 
pinning regime for a certain temperature range, above 
which the peak effect and the onset of plasticity simul- 
taneously appear. In our model the peak effect is caused 
not by the onset of plasticity in the a — b plane but by 
the onset of plasticity in the c-axis [fi"7f with the plasticity 
only occurring for T > Td c - We find that the tempera- 
ture Tdc moves further below T m as the disorder in the 
sample is increased. Thus for large enough disorder the 
2D phase extends all the way to T — 0, and the dynamic 
phase diagram becomes the same as the one proposed in 
in Ref. f§. 

In Fig. 4 we show the V x versus fd curves, which cor- 
respond to experimental V(I) curves, for increasing T. 
For T < Tdc there is a sharp elastic depinning transition 
at which all the vortices start moving at once in the 3D 
state. For T > Td c the depinning response is more con- 
tinuous since the 2D depinning occurs inhomogeneously, 
with certain regions moving while other regions remain 
pinned. As T is further increased in the 2D regime the de- 
pinning transition falls at driving forces lower than that 
at which the 3D depinning transition occurred. Although 
for these higher temperatures the 2D pancakes depin at 
a lower applied drive than the 3D state, the overall veloc- 
ity above depinning is still less than that of a coherently 
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moving 3D system at the same drive so that a crossing 
of the IV curves occurs. 

To summarize, we have numerically investigated the 
dynamics of magnetically coupled vortices in layered su- 
perconductors interacting with quenched disorder. For 
sufficiently strong disorder the depinning is 2D. Here the 
vortices flow plastically in each layer and arc uncorrclatcd 
from layer to layer. For higher drives there is a sharp 
hysteretic transition to a coherently moving ordered 3D 
vortex lattice. For weak disorder there is a non-hysteretic 
depinning transition from a pinned 3D to a moving 3D 
state. The two stage 2D-3D hysteretic depinning tran- 
sition is in good agreement with the theoretical predic- 
tions of Ref . Jll[] . We have also mapped out the dynamic 
phase diagram as a function of driving force and temper- 
ature. At low temperatures we observe non-hysteretic 3D 
pinned to 3D moving transitions. For increased temper- 
ature there is a static 3D-2D transition which coincides 
with a peak in the critical current as well as the onset of 
plasticity in the c-axis. The drive at which a dynamically 
induced 2D-3D dynamic transition diverges as T m is ap- 
proached. Our dynamic phase diagram is very similar to 
that proposed in Ref 0. 
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FIG. 3. The fd versus T/T m dynamic phase diagram for 
a system with f p — 0.01/q. Filled circles: Depinning line. 
Open squares: Dynamic reordering line. 
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FIG. 4. V x versus fd responses for varied T. Cir- 
cles: temperatures below the decoupling temperature Td c : 
T/T m = 0.00027 and 0.0013. Lines, right to left: T/T m = 
0.005, 0.014, 0.05, 0.11, 0.22, 0.5, 0.9. 
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